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Abstract 

It is known that the so-called Bercovici-Pata bijection can be explained in terms of certain random 
matrix ensembles {Md) d>1 whose asymptotic spectral distributions are free infinitely divisible. We in- 
vestigate matrix Levy processes with jumps of rank one associated to these random matrix ensembles 
introduced in [5| and [8|. A sample path approximation by covariation processes for these matrix Levy 
C*") ■ processes is obtained. As a general result we prove that any d x d matrix subordinator with jumps of rank 

one is the quadratic variation of an Revalued Levy process. In particular, we have the corresponding 
result for the matrix subordinators with jumps of rank one associated to the random matrix ensembles 
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AMS 2010 Subject Classification: 60B20; 60E07; 60G51; 60G57. 

^ ; 1 Introduction 
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New models of infinitely divisible random matrices have emerged in recent years from both applications and 
theory. On the one hand, they have been important in multivariate financial Levy modelling where stochastic 
volatility models have been proposed using Levy and Ornstein-Uhlenbeck matrix valued processes; see [2J, 
[3], [1] and |12j . A key role in these models is played by the matrix subordinators and more general matrix 
covariation processes. 

On the other hand, in the context of free probability, Bercovici and Pata [7] introduced a bijection A from 
the set of classical infinitely divisible distributions to the set of free infinitely divisible distributions. This 
bijection was explained in terms of random matrix ensembles by Benaych- Georges [5] and Cabanal-Duvillard 
[5] , providing in a more palpable way the bijection A and producing a new kind of infinitely divisible random 
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matrix ensembles. Moreover, the results in [5] and [5] constitute a generalization of Wigner's result for the 
Gaussian unitary ensemble and give an alternative simple infinitely divisible random matrix model for the 
Marchenko-Pastur distribution, for which the Wishart and other empirical covariance matrix ensembles are 
not infinitely divisible. 

More specifically, it is shown in [5 and [8 that for any one-dimensional infinitely divisible distribution /i 
there is an ensemble of Hermitian random matrices (Af<j)d>i, whose empirical spectral distribution converges 
weakly almost surely to A(/n) as d goes to infinity. Moreover, for each d > I , Md has a unitary invariant 
matrix distribution which is also infinitely divisible in the matrix sense. From now on we call these models 
BGCD matrix ensembles. It is possible to obtain BGCD models of symmetric random matrices rather than 
Hermitian. Indeed, slight changes in the proof of [3 Theorem 3.1] give for each d > 1, a d x d real symmetric 
random matrix Md with orthogonal invariant infinitely divisible matrix distribution. We consider additional 
facts of BGCD models in Section 3. 

A problem of further interest is to understand the matrix Levy process {Md(t)} t>0 associated to the 
BGCD matrix ensembles. It was pointed out in [5], [IT] that the Levy measures of these models are 
concentrated on rank-one matrices. This means that the random matrix Md is a realization, at time one, of 
a matrix valued Levy process {Md(t)} t>0 with rank one jumps AMg(t) — Md(t) — Md(t~). 

The purpose of this paper is to study the structure of a d x d Hermitian Levy process {Ld(t)} t>0 
with rank-one jumps. It is shown in Section 4 that if Ld is a real symmetric nonnegative definite matrix 
Levy process, it is the quadratic variation of an Revalued Levy process Xd, being the converse part of 
the folklore. The process Xd is constructed via its Levy-Ito decomposition. In Section 5 we consider new 
realizations in terms of covariation of C d -valued Levy process for matrix compound Poisson process as well 
as sample path approximations for Levy processes associated to general BGCD ensembles. A new insight 
on Marchenko-Pastur's type results for empirical covariance matrix ensembles was recently given in [5] by 
considering compound Poisson models (then infinitely divisible). In this direction our results show the role 
of covariation of d-dimensional Levy process as an alternative to empirical covariance processes. 

For convenience of the reader, and since the material and notation in the literature is disperse and 
incomplete, we include Section 2 with a review on preliminaries on complex matrix semimartingales and 
matrix valued Levy processes that are used later on in this paper. 

2 Preliminaries on Matrix Levy processes 

Let Mdxq = Mdxij (C) (respectively M^xg (R)) denote the linear space of dx q matrices with complex (respec- 

1/2 

tively real) entries with scalar product (A, B) — tr (AB* ) and the Frobenius norm || A|| = [tr (AA*)Y< where 

1/2 

tr denotes the (non normalized) trace and |A| = (AA*) ' . If q = d, we write = M^xd- The set of Her- 
mitian random matrices in is denoted by H^. Likewise, let U^xg = U^xg (C) = {U € M.d Xq '■ U*U = I q } . 
If q = d, Vd = Udxd- We denote by Hjm the set of matrices in of rank I and by (H^~) the set of 
positive (nonnegative) definite matrices in Md- Likewise = M d (i) n H^. 

Covariation of complex matrix semimartingales An M pxg -valued process X = 
{(xij)t} t>0 is a matrix semimartingale if (xij)t is a complex semimartingale for each i = 1, ...,p, j — I, q. 
Let X — {(xij)t} t>0 £ M pxq and Y — {{yij)t} t>0 € M qxr be semimartingales. Similar to the case of 
matrices with real entries in [2], we define the matrix covariation of X and Y as the M pxr - valued process 
[X,Y] := ([X,Y] t ) t>0 with entries 




(1) 
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where [xik,ykj] t is the covariation of the C-valued semimartingales {(xik)t} and {(xkj)t}l see [131 pp 83]. 
One has the decomposition into a continuous part and a pure jump part as follows 

[X, Y] t = X Y + [X c , Y\ + J2 (AX.) (AY,) , (2) 

s<t 

where [X c , Y%° := Yl =1 [x lk ,y ko ] c t . 

We will use the facts that [X] — [X, X*} is a nonnegative definite pxp matrix, that [X, Y] T = [Y T , X T1 \ 
and that for any nonrandom matrices A g M mxp , C € M rxn and semimartingales X S M. pxq ,Y £ M qxr 

[AX, YC] = A [X, Y] C. (3) 

The natural example of a continuous semimartingale is the standard complex d x q matrix Brown- 
ian motion B — {B t } t>Q = {{bji(t))} t>0 consisting of independent C-valued Brownian motions bji(t) = 
Ke(bji(t)) +ilm(bji(t)) where Ke(bji(t)), Im(bji(t)) are independent one-dimensional Brownian motions with 
common variance t/2. Then we have [B,B*]l J = J2l=i [bik,bjk] t = qtSij and hence the matrix quadratic 
variation of B is given by the d x d matrix process 

[B,B\ = qtl d . (4) 

The case q = 1 corresponds to the C d -valued standard Brownian motion B. We observe this corresponds to 
[B, B*]. = tld instead of the common 2tl d used in the literature. 

Other examples of complex matrix semimartingales are Levy processes considered next. 

Complex matrix Levy processes An infinitely divisible random matrix M in M^xg is characterized by 
the Levy-Khintchine representation of its Fourier transform Ee ltr( - e M ' = cxp(?/>(0)) with Laplace exponent 

V>(6) = itr(6** ) - \tv (e*AG*) + f I e itr ( e *«) - 1 - ) v(dO, 6 e M dxq , (5) 

1 JM dxq \ 1 + HCII / 

where A : M ?X( i — > Mdxq is a positive symmetric linear operator (i.e. tr (Q* A$*) > for $ g Mdxq and 
tr (Q*.A&1) = tr (QIAQ*,) f° r ©i,©2 G M d x q ), v is a measure on M dxq (the Levy measure) satisfying 
^({0}) = an d / M (INI A l)z^(dx) < oo, and ^ € Md X? - The triplet (.A, ^, vp) uniquely determines the 
distribution of M. 

Remark 1 TTie notation AO* means the linear operator A from M qxd to Mdxq acting on O* € M qX d- 
Some interesting examples of A and its corresponding matrix Gaussian distributions are: 

(a) AO* — O. This corresponds to a Gaussian matrix distribution invariant under right and left unitary 
transformations in\Jd andV q , respectively. 

(b) AO* = E1OE2 for Ei G and £2 G M„ . In this case the corresponding matrix Gaussian 
distribution is denoted by Ndxg(0, £1 <g> £2) and £1 (g> £2 is called a Kronecker covariance. It holds that if 
N has the distribution Nrf X9 (0,I,j ® \ q ), then T, 1 ^ 2 NT,^ 2 has distribution Nd Xg (0, £1 (g> £2). 

(c) When q = d, AO* = tr(9)Irf is the covariance operator of the Gaussian random matrix gld where g 
is a one- dimensional random variable with a standard Gaussian distribution. 

Let E>dxq the unit sphere of Md X? and let M2 xg = Md Xg \{0}. If v is a Levy measure on Sdxq, then there 
are a measure A on Sdxq with \(§dxq) > and a measure for each £ € Sdxq with ^((0, 00)) > such that 

v(E) = [ A(d£) f l £ «)^(du), E £ B(M° dxq ). 

JS dxq J(0,oa) 

We call (A, v^) a polar decomposition of v. When d = q = 1, ^ is a Levy measure on R and A is a measure 
in the unit sphere §i x i = { — 1, 1} of K.. 
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Any M^xq- valued Levy process L = {Lt} t>0 with triplet (A, is, is a semimartingale with the Levy-Ito 
decomposition 



L t = m + B A (t) 



[0,t] 



xJl{&s, dx) 



il<i 



[0,t] «/||x||>l 



a;Ji(ds, dx), £ > 0, 



(G) 



where: 

(a) {B A (t)} t>Q is a Mrf X(? -valued Brownian motion with covariance A, i.e. it is a Levy process with 
continuous sample paths (a.s.) and each B A (t) is centered Gaussian with 



E {tr(6*-B^(i))tr (0*.B A (s)) } = min(s, t)tv (0\AO* 2 ) for each 9i,e 2 G 



(b) Jl(-) is the Poisson random measure of jumps on [0, oo)xM° That is, Jl(£, E) = # {(0 < s < t : AL S £ E} , 



E e 



with intensity measure Leb® v, and independent of {B A {t)} t>0 , 



(c) Jjj is the compensator measure of Jl, i.e. 

Ji(di,dx) = JL(dt,dx) — dtv(dx); 

see for example pQ for the most general case of Levy processes with values in infinite dimensional Banach 
spaces. 

An Mrf X? - valued Levy process L = {L t } t>0 has finite variation if and only if its Levy-Ito decomposition 
takes the form 

L t = m + I I xj L (ds,dx) = m + Val s , t > o, (7) 

where $0 = ^ — i||a:l|<i xv(dx). 

The matrix quadratic variation ([2} of L is given by the Mj"-process 



[L] t = [B A ,B* A ] t + [ f xx*J L (ds,dx) = [B A ,B A ] t + y2AL s AL* s . 
In Section 3 we prove a partial converse of the last result in the case q = 1. 



(8) 



Remark 2 On the lines of Remark]]] we have the following observations for the quadratic variation of the 
continuous part in @): 

(a) When AO* = O, [B A ,B* A ] t = gild- This follows from |^J) since B A (t) is a standard complex d x q 
matrix Brownian motion. 



(b) When AO* = Ei9E 2 for £j € M f j and E 2 6 



ffeaf £4(4) = Ej /2 B(t)E2 /z w/iere 



a/2 



_B = {Bt} t>0 is a standard complex d x q matrix Brownian motion. Then, using $3ty we /iaue 



.4Ji 



E^E^E^E^ 2 



-1I/2 



BE^E^B* 



1I/2 



itr(E 2 )S 1; 



where we have also used the easily checked fact 



beI /2 ,eI /2 b* 



*tr(E 2 ). 



fcj When q = d and AO* = tr(0)ld, we have [B A ,B A ] t = tld since Bj,(t) = b{t)Id where b = {b(t)} t>0 
is a one- dimensional Brownian motion. 



Approximation of Levy processes The following are useful results on the sample path approximation 
of complex matrix Levy processes; see fTUl Th 15.17] and [TH Th. 8.7]. They follow from their corresponding 
real vector case by the usual identification of Md Xg — > M. 2dq via A — > vec(A), A £ M.dxq and the fact that 
tr (A*B) = vec(A)*vec(B), where vec(A) is the dq column complex vector obtained by stacking the columns 
of A one down each other. 
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Proposition 3 Let L,L ,L , ... be complex matrix Levy processes in M^xg with L n (l) A L(l). Then th ere 
exist processes L n with the same distribution that L n such that 



sup 

0<s<t 



VI - L< 



Pr 



Vt > 0. 



Proposition 4 Let M n , n = 1,2, ... be infinitely divisible random matrices in Mdxq with triplet (A n , v n , Vf"). 

Lei M be a random matrix in Wldxq- Then M n — » M i/ and on/y i/ M is infinitely divisible whose triplet 
(A, v 1 satisfies the following three conditions: 

a ) If f '■ Mdxg — > ^-dxq is bounded and continuous function vanishing in a neighborhood of then 



lira 

n— >oo 



/(£K(d£) = / /(0"(d0 



&J Define the positive symmetric operator A n ' e : M. qxc i — > M^xg by 

ti(e*A n - e e*) = tr(e*A n e*)+ [ |tr(e*£)| 2 i/ n (d0 /< 



'or € 



limlimsup|tr(eM n ' e e*)-tr(eMe*)| =0, /or6eM d 



eJ.0 



3 BGCD random matrix ensembles 

We now consider the matrix Levy processes associated to the BGCD matrix ensembles {Md)d>i mentioned 
in the introduction. 

When [i is the standard Gaussian distribution, M4 is a Gaussian unitary invariant random matrix, A(/j,) 
is the semicircle distribution and {Md(t)} t>0 is the Hermitian matrix-valued process given by Md{t) = 
(l/vd+l) (Bt + dg t ld) where B t is a d x d Hermitian matrix Brownian motion independent of the one- 
dimensional Brownian motion g t ; see Remark 3.5]. 

Likewise, if \x is the Poisson distribution with parameter A > 0, {Md(t)} t>Q is the dxd matrix compound 

Poisson process Md(t) — Ylk=i u k u k* where { u t}k>i * s a se< t uence of independent uniformly distributed 
random vectors on the unit sphere of C d independent of the Poisson process {N t } t>0 , and A(/i) is the 
Marchenko-Pastur distribution of parameter A > 0; see [5J Remark 3.2]. We observe that in this case 
{Md(t)} t>Q is a matrix covariation (quadratic) process rather than a covariance matrix process as in the 
Wishart or other empirical covariance processes. 

Proposition [5] below collects computations in [5] , [8] and [9] to summarize the Levy triplet of a general 
BGCD matrix ensemble in an explicit manner. We write §(HLi(x)) = HWi) n §dxd- 

Proposition 5 Let p. be an infinitely divisible distribution in K with Levy triplet (a 2 ,ip,is) and let (Md)d>i 
be a BGCD matrix ensemble for A(/x). Then, for each d > 1 Md has the Levy-Khintchine representation (0) 
with Levy triplet (Ad, ^d, v<i) where 

a) * d = iPdld 

b) 

^e = a 2 -i-(e + tr(e)I d ),eeH d . (9a) 
a + 1 

c) 

v d {E)=d( [°° l E (rV)v v (dr)n(dV), E € B (H d \ {0}) , (10) 

JS(M d(1) ) JO 
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where vy = v\(n, oo) or H(-oo,o) according to V > or V < and II is a measure on S(M d ^) such that 

11(D) = J J l D (tV)\(dt)u d (dV), D e B (S(U d{1) )) , 
S(H d(1) )nl+ {-I. 1 } 

where A is the spherical measure of v and tu d is the probability measure on S(Hd(i)) f~l given by the 
probability measure on the set of matrices of rank-one induced by the transformation u — > V = uu* , where u 
is a uniformly distributed column random vector in the unit sphere of C d . 

Proof, (a) It follows from the first term in the Levy exponent of M d in page 635 of [8], where the notation 
A d (n) is used for the distribution of M d . For (b), the form of the covariance operator A d was implicitly 
computed in the first example in Section II. C of [8]. Finally the polar decomposition (|TU1) was found in [5]. 
■ 

The Levy-Ito decomposition of the Levy process associated to the BGCD model M d is given by 



M d (t) = ijjtdl d + B A (t) + / xJ d (ds,dx) + / / xJ d (ds,dx), (11) 

J lo,t] ^{||x||<i}nH d(1) J[o,t] J {\\x\\>i}rM dll) 

where * > 0, A d Q = a 2 ^{Q + tr(6)I d ), J d (t,E) = #{0<s<t: AM d (s) e E} = J d (t, E n H d(1) ) for any 
measurable E G MK, Its quadratic variation is the matrix subordinator 

[M d ] t = a 2 tl d + [ [ xx*J d (ds, dx) = a 2 tl d + V AM d (s) (AM d (s))* . 
•V*] Jm° w s < t 

4 Bounded variation case 

It is well know that the quadratic variation of a one-dimensional Levy process is a subordinator. The 
following result gives a converse and a generalization to matrix subordinators with rank-one jumps. The one 
dimensional case is given in [15( Lemma 6.5]. 

Theorem 6 Let L d = {L d (t) : t > 0} be a Levy process in M. d (R) whose jumps are of rank one almost 
surely. Then there exists a Levy process X = {X(t) : t > 0} in M d such that L d (t) = [X] t . 

Proof. We construct X as a Levy-Ito decomposition realization. For each d > 1, L d = {L d (t) : t > 0} is an 
(R)-process of bounded variation with Levy-Ito decomposition 



L d (t) = ^ t+ / xJ Ld (ds,dx),t>0, 
J[0,t] Jm+ (1) (r)\{o} 

where \I/o G ^- d (^) an d Jh d is the Poisson random measure of L t . Let Vl 4 be the intensity measure of Ji, d . 

Let Cf — {x — (xi,x<z, . . . , x n ) : x\ > 0} and C| = {x = (xi, x 2 , ■ ■ ■ , x n ) : X\ < 0}. Note that Cf and C| 
form a partition of R d . Let (fj : M + x M+ (1) (R) -> M + x be defined as (i, V) = (i, x) where V = xx T 

and x E Cf for j = 1, 2. Let ^ : M+ (1) (R) -> be defined by ^ (V") = x for V" = xx T and x £ Cf for 
3 = 1,2. 

Let us define (ds, dx) = ( Ji d o ^ ) (ds, dx) for j = 1,2, the random measure induced by the transfor- 
mation ifj which is a Poisson random measure on R + x Cf. Observe that E \J 3 x (t, A)\ = E j£ d o ^t 1 ({*} x A) 

tvi, d (jpj 1 (A)) — t {yi, d ° Vj 1 ) {A) for A <E B (Cj\ {0}). Let us denote v x = VL d o^pj 1 which is a Levy mea- 
sure on Cf since 



(l A |x| 2 ) u x (dx) = / (l A |x| 2 ) ^ o ^(dx) = / (1 A tr (rf)) „ £- o ^(ds) 
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(lAtr(V))(v Ld o^)ori(dV) 



+ 

d(i) 



(lAtr(V))i/ £(l (dK) <oo, 



where \yh d ° ^ ) ° f^ 1 = VL d , with / (x) = xx 1 and we have used that tr (V) < a ||V|| for some constant 
a where ||V|| denotes the Frobenius norm of V. Thus v 3 x is the intensity measure of the Poisson random 
measure Jl . 



Let us define Jx(ds,dx) — (hJx + h^x) (ds, da;) which is a Poisson random measure on 



Now 



let us define the measure vx — \ v L d ° Pi 1 + \ v L d ° ¥2 ■ Observe that vx is a Levy measure on R d which 
satisfies E [Jx(t,A)} — tvx {A) for A G B (R d \ {0}). Hence fx is the intensity measure of Jx- 
Let us take the Levy process in R d 



X(t) = \* \ 1/2 B A (t) 



[o,t] JR d n{|x|<i} 



xJx(ds, dx) 



[o,t] Js. d n{\x\>i} 



xJ x (ds,dx),t>0, (13) 



B A is a Revalued standard Brownian motion with covariance A. Then its quadratic variation is given by 



[X] t =[\^ 2 B A ,B* A \^ 2 
= Vatl d - 



xx*Jx(ds, dx) 



[0,t] JR d \{0} 

' I 

[o,t] Jm d \{ 



[0,t] JR d \{0} 

xx 1 (\J X + h J x) (ds,dx) 



xxrj Ld o^ {ds , d x) + U [ xxrj Ld o^ {ds , d x) 

0} J[0,t] JR d \{0} 



= *o< + I / / V ( J Ld o (^r 1 ) h- 1 (ds, dV) 

J[o,t] Jm+ (1) \{o) 

+ §/ / ViJL.o^oh-^dV) 
J[o,t] Jm+ (1) \{o} 

= *o* + | / / FJ i(J (ds,dF) + ± / / V/J Ld (ds,d^) = L t , 

J[o,t] -/m+ (1) \{o} i[o,t] Jm+ (1) \{o} 

where (JL d fj 1 ) ^ 1 — ^L d , with ft, (t, x) = (t, xirJ) . ■ 
For the general bounded variation case we have 

Theorem 7 Let Ld = {Ld(t) : t > 0} be a Levy process in M.d (R) of bounded variation whose jumps are of 
rank one almost surely. Then there exist Levy processes X — {X(t) : t > 0} and Y = {Y(t) : t > 0} in R d 
such that L d {t) = [X] t - [Y] t . 

Proof. For each d > 1, Ld = {Ld{t) : t > 0} is an MLj (R)-process of bounded variation with Levy-Ito 
decomposition 

L d (t) = m+ I I xJ Ld (ds,dx),t>0, (14) 

J[0,t] J%)(8)\{D} 

where * G M d (R) and JL d is the Poisson random measure of L t . Observe that L d = L d — L 2 d where L\ and 
L d are increasing Levy processes in M d (R), if we take their corresponding Levy-Ito decompositions with 
*i > 0, J\ d and f 2 > 0, J 2 Ld then * = #1 - * 2 - 

Following the terminology in the proof of Theorem [S] we construct Levy processes X and Y in R d as 
Levy-Ito decompositions with same covariance A, 

X(t) = \V 1 \ 1/2 B A {t)+ [ [ xJ x {ds,dx)+ [ [ xJ x {ds 7 dx),t>0, 

J[0,t] JR d n{\x\<l} J[0,t] JR d n{\x\>l} 
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Y(t) = \V 2 \ 1/2 B A (t) + [ [ xJ Y (ds,dx)+ [ [ xJ Y (ds,dx),t>0, 

J[0,t] JR d n{\x\<l} J[0,t] JR d n{|a;|>l} 

where Jx(ds,dx) = (\J X + \ J~x) (ds, dx) is constructed in terms of tpj ■ : R x MTt^ — > C d and the Poisson 
random measure J Y (ds,dx) on R d is constructed similarly by using the induced measure JL d (t, ■) ° Pj (■) 
for the transformation /3j : E X — > C d defined as f3j (t, V) = (t,x) where V — —xx T and x £ C d for 
j = 1,2. 

Observe that 

[X] t = #!*+/ / ra T J x (ds,da;) = *ii+ / / yj Lti (ds,d^) 
7[o,t] Vr«\{o} 7[o,t] Jm+ (1) \{o} 

and 

[Y] t = -$f 2 t+ / xx T J Y (ds,dx) = y 2 t- / (-.to 1 ) Jy(i,dar) 

J[0,t] JR d \{0} J[0,t] JC d \{0} 

= ¥2*- / / 7J id (ds,dV0. 

J[0,t] JMJ (1) \{0} 

The result follows from (fT4|) . ■ 

Next we consider the matrix Levy processes associated to the real BGCD matrix ensembles (Md)d>i 
mentioned in the Introduction. In this case for each d > 1, Md is a d x d real symmetric random matrix 
whose distribution is invariant under orthogonal transformations and infinitely divisible in (R) with Levy 
measure concentrated in rank one matrices. 

We have the following two consequences of the former results. 

Corollary 8 Let Md = {Md(t) : t > 0} be the matrix Levy process associated to the real BGCD random 
matrix ensembles. 

a) Let /i be the infinitely divisible distribution with triplet (0, ip, v) associated to Md such that f\ x \ <1 (1 A x) v (da;) < 
00, 00, 0]) = and ipo :— ip — J x<1 xv{dx) > 0. Let us consider the Levy-Ito decomposition of Md(t) in 



M d (t) = tpotd + / xJ Md {ds,dx). 

J[0,t] JM+ (1) (R)\{0} 

Then there exists a Levy process X = {X(t) : t > 0} in R d such that M d {t) = [X] t , where 

X(t) = \ifj \ 1/2 B A (t)+ I f xJ x {ds,dx)+ I f xJ x (ds,dx), t > 0, 

J[0,t] JC A n\{\x\<l} J[0,t] JC d n\{\x\>l} 

Ba is a W 1 -valued standard Brownian motion with covariance A, the Poisson random measure Jx is given 
by J x {ds,dx) = (\JM d ° ^ + \Ju d ° f^ 1 ) {ds,dx). 

b) If Md has bounded variation then there exist Levy processes X = {X(t) : t > 0} and Y = {Y(t) : t > 0} 
in R d such that M d (t) = [X] t - [Y] t . 

Remark 9 The idea of the proof of the above results does not apply to the case of Hermitian random 
matrices. This is due to the fact that the unit sphere in C is uncountable, while in R it has two points. More 
specifically, a decomposition of C d in a finite or countable union of sets C d as was done for M. d in the proof 
of Theorem [S| is not longer possible. The following section presents some other covariation representations 
for the Hermitian case. 
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5 Covariation matrix processes approximation 

We now consider approximation of general BGCD ensembles by BGCD matrix compound Poisson processes 
which are covariation of C d -valued Levy processes. 

The following results gives realizations of BGCD ensembles of compound Poisson type as the covariation 
of two C d -valued Levy processes. Its proof is straightforward. 

Proposition 10 Let fi be a compound Poisson distribution on R Levy measure v and drift tp € R and let 
(Md)d>i be the BGCD matrix ensemble for A(/i). For each d > 1, assume that 

i) ((3j)j>i is a sequence of i.i.d. random variables with distribution ^/^(R). 

ii) {uj)j>i is a sequence of i.i.d. random vectors with uniform distribution on the unit sphere ofC d . 
Hi) {N t } t>Q is a Poisson process with parameter one. 

Assume that (Pj)j>i, (uj)j>i and {N t } t>0 are independent. Then 

a) Md has the same distribution as Md(T) where 

N t 

M d (t) = ipi d +j2 1 ^ °- ( 15 ) 

b) M d (-) = [X d ,Y d ]{.) where X d = {X d (t)} t > , Y d = {Y d {t)} t > are the C d -valued Levy processes 

N t 

X d (t) = VW\Bt + ^2yJ\l3j\uj, t>0, (16) 

Y d {t) = sign (V) VW\Bt + si § n (ft") >J\Pi\ u i> 1 ^ °. ( 17 ) 

and B — {B t } t>a is a C d -valued standard Brownian motion independent of (Pj)j>i, (%)j>i and {N t } t>Q . 

For the general case we have the following sample path approximation by covariation processes for Levy 
processes generated by the BGCD matrix ensembles. 

Theorem 11 Let \i be an infinitely divisible distribution on R with triplet (a 2 ,-0,^) and let (Md)d>i be the 
corresponding BGCD matrix ensemble for A(p). Let d > 1 fixed and assume that for n > 1 

i) (/3™)j>i is a sequence of i.i.d. random variables with distribution 

ii) (w")j>i is a sequence of i.i.d. random vectors with uniform distribution on the unit sphere of C d . 
Hi) N n — {-/V"} t>0 is a Poisson process with parameter n. 

iv) B n = {B"} t>0 is a C d -valued standard Brownian motion. 

v) (/3")j>i, (Uj)j>i,N n and B n are independent. 
Let 

A? 

X2{t) = y/W\Bt+Y,V\l 3 j\ U i> *^ ' ( 18 ) 
j=l 

AT 

Y?(t) = sign (V) y/W\&? + E si S n (^") \fW\ u l' 1 ^ °- ( 19 ) 

3 = 1 

Then for each d > 1 there exist M d (C) -valued processes M d l = (m£(*)) such that M™ = [X%, Y 7 d l ], 

> d>l 

M2(s) - M d {s) 

where {Md(t) : t > 0} is the M<j (C)-valued Levy process associated to (Md)d>i- 



sup 

0<s<t 



A 0, Vi > 7 
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Proof. By the compound Poisson approximation for infinitely divisible distributions on K (see for example 
[141 pp 45]), let n n be an infinitely divisible distribution such that fi n — > fi, where we take the triplet of y, n 
as (A n ,^ n ,v n ), A n = 0, V" = J i+b* 1 '" ( dx ) and yn = n t i *"- 
Following the proof of Theorem [TOl on the same lines we have that 



M":= [X2,Yr) t = iPtI d + J2P 



n u n u n * 

3 3 3 ' 



which is a matrix value compound Poisson process with drift whose triplet (A2,ip^,v^) is given by A\\ = 
0, = and 

v2{E)=df H l E (rV)v v (dr)U(dV), E € B (W d \ {0}) , (20) 

JS(M da) ) JO 



where v v and II are as in JI 

We will prove that M% — >■ M d by showing that the triplet (A^, V$ , v2) converges to the triplet (Ad, i[>d, v d) 
of the BGCD matrix ensemble in Proposition [5] in the sense of Proposition [4] 
We get ip^Id = ipld — > ipld as n — > oo. 

Let / : Hrf^) — >■ 1 be a continuous bounded function vanishing in a neighborhood of zero. Using the polar 
decomposition (fT0|) for we have 



/ (0 v2 <M) =d f (rV) v v (dr) U (dV) 

Hd(i) Js(M d(1) )J0 



= d / f(trV)v v (dr)\ n (dt)w d (dV). (21) 

JS(H d(1) )nH+ 7{-i,i}Jo 



For V € S(H d(1) ) n Hrf fixed, 

/ / f(trV)v v (dr)\ n (dt) = \ n ({l}) f(rV)v n (dr) 
J {-1,1} Jo Jo 

+ A"({-1}) f° f{rV)v n (dr) 



and as a function of r, / (rV) is a real valued continuous bounded function vanishing in a neighborhood of 
zero, for which 



A" ({1}) / / (rV) P n (dr) — » A ({1}) / / (rV) v (dr) 
Jo Jo 

and 

A" ({-1}) / / (r7) i/ B (dr) — ► A ({-1}) / / (rF) i/ (dr) 

J — oo ./ — OO 

Then from (HD 



/ (0 "2 (dt) -^d f (trV) v v (dr) A (dt) u d (dV) 

H d( i, JS(H d(1) )nHJ J{-1,1} Jo 

/(rV)^(dr)n(dl/) = / f(i)v d (dO- 



Next, we verify the convergence of the Gaussian part. 



10 



Let us define, for each e > and n > 1, the operator A n ' £ : — > M d by 

tr(OA"> s e) = [ \tr(QO\ 2 Vd(dO- 

JH\\<e 

From (H| we get 

/ |tr(90| 2 ^(^) 

JllfIKe 



rf / / l{||ry||<e}(^)|tr(rey)| 2 ^(dr)n(dy) 

Js(U d(1) ) Jo 



oo 

2 1+^ ^m|2 . 



d/ / / l {r < e }{rtV)r 2 \ti{QV)\ z ^{dr)\{dt)uj d {dV) 



ld ( i))nH+ J {-1,1} Jo 

2 



d _ / l{r<e}(rV)r z \tr(QV)\ z v n (dr)Lu d (dV) 

JS(H d(1) )nH+ 



= d/ _ |tr(6V0r / r z v n {dr)uj d {dV). 

Since Ji r i <E . r 2 v n (dr) — > a 2 then 

|tr (60I 2 v n d (d£) — -> da 2 ^„ |tr (9im*)| 2 

^(tr(e 2 ) + (tr(9)) 2 ), 



ll?ll<£ 

da 2 



which gives the operator in part (b) of Proposition [5] Thus M d l — ► M d . Finally the conclusion follows from 
Proposition [3l ■ 

Acknowledgement 12 This work was done while Victor Perez- Abreu was visiting Universidad Autonoma 
de Sinaloa in January and May of 2012. 
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